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GGSS  RAIL  GARRISON  AND  VAN  SURVEY  EXPERIENCE 


Operational  problems  are  described  with 
observations  and  results  relating  to  the 
acceleration  environment  on  the  rail  car. 
Techniques  used  to  identify  and  compensate 
for  sensitivities  are  covered.  Survey 
results  are  presented. 
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ABSTRACT 


By 

ANDREW  D.  GRIERSON 

BELL  AEROSPACE  TEXTRON 
Division  of  Textron,  Inc. 
P.O.  Box  One 

Buffalo,  New  York  14240-0001 


GRAVITY  GRADIOMETER  SURVEY  SYSTEM 
STAGE  I  DATA  REDUCTION 


Suppression  of  motion- induced  signals  and  errors  is  vital 
to  the  measurement  of  gravity  gradient  signals  on  a  moving 
vehicle.  Post-mission  processing  of  data  from  the  Bell 
Gravity  Gradiometer  Survey  System  (GGSS)  uses  intrinsic 
features  of  the  system  design  to  accomplish  this. 

GGSS  mounted  on  a  railroad  freight  car  experienced  a 
high  vibration  environment.  Post-mission  compensation 
of  the  recorded  data  has  removed  much  of  the  sensitivity 
to  this  environment. 

Improvements  in  platform  control  resulted  from  experiences 
in  railroad  operation.  Subsequent  road  trials  have  demon¬ 
strated  the  ability  to  measure  gradients  of  local  topographic 
features. 
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MODIFIED  PROCESSING  FOR  RAIL  SURVEY 
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RMS  ACCELERATION  SENSITIVITY 
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GRADIENTS  MEASURED  ON  A  SECTION  OF  ROAD  SURVEY  SHOWING  LOCAL  TOPOGRAPHIC  FEATURES 
5  PASSES  PERFORMED  OVER  A  PERIOD  OF  2  DAYS. 


WHAT  HAS  THE  EXPERIENCE  OF  OPERATION  OF  THE  GGSS  IN  THE  AIR  AND  ON 
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THE  AGE  OF  MOVING  BASE  GRAVITY  GRADIFNT  MFAf^l 


iU; 


a:: 

o 

u. 


=31 

cs 

00 

<: 

zc 


=3  LJ 
00 

UJ  UJ 
QC  O 


O 

ac 

o 


r-  os 


a 


00 

UJ  “ 

g;  ^ 


(A 


> 


> 

c 

as 

09 

09 


00 


oe 

09 

UJ 

< 

>“ 

09 

OJ 

•-M 

«c 

> 

• 

cc 

u_ 

< 

00 

3 

o 

UJ 

OJ 

OJ 

09 

UJ 

a 

< 

z 

> 

9Z 

o. 

00 

zc 

09 

Q- 

00 

OJ 

=C 

2 

UJ 

1— 

36 


ABSTRACT 


By 

ALBERT  JIRCITANO 

BELL  AEROSPACE  TEXTRON 
Division  of  Textron,  Inc. 

P.O.  Box  One 

Buffalo,  New  York  14240-0001 

GRAVITY  BASED  PASSIVE  NAVIGATION 


Passive  covert  Inertial  Navigation  System  (INS) 
updating  can  be  implemented  by  comparing  measured 
gravity  gradients  with  mapped  values  and  using  the 
error  in  an  optimal  filter  to  define  in  real  time 
the  INS.  A  parametric  study  is  carried  out 
showing  performance  as  a  function  of; 

0  Gradiometer  accuracy  and  number  of  gravity 
gradient  sensors. 

0  INS  -  gyro  and  accelerometer  accuracy 
0  Map  quality 

0  Gravity  field  characterization  (field  intensity 
and  frequency  content) 

0  Altitude 

0  Velocity 
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GRAVITY  GRADIENT  PASSIVE  NAVIGATION 
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CAVITY  BASED  PASSIVE  NAVIGATION  CHARACTERISTICS 


o 

g  a 

UJ  oc 


<  ca£ 
CO  o 

>  OS 

<  LU 


rs  c<o 
o  o 

.«5  Q„ 


CO 

LU 

oe 

«E 

53 

ce: 

C7 

LU 

LU 

Qi 

< 

a. 

* 

g 

OS 

UJ 

o 

> 

o 

LU 

LU 

< 

i"< 

QC 

:> 

LU 

< 

o. 

os 

o 

CO 

• 

• 

CO 

OS 

* 

^  g 
I  hi 


41 


EXTENSIVE  DMA  AND  NOO  DATA  BASES. 


Summary  of  Parametric  Study  Conditions 


U 

Ui 

*/> 

Q 

< 

oe 

LU 

o 

o 

o 

o 

o 


I 

UJ 

«/l 

5 

Z 

UJ 

H 

X 


UJ 

H* 

I/) 

> 


i/t 

UJ 

Q 

O 

s 


o 

o 


o 

o 

<N 


S  ^ 

“  <a 

Ui  m 

S  ^ 

O  Z 

5  < 

<  !" 
oe  Z 

V9  O 


o 
o 
Q  ^ 

ui 

“■  UJ 

^  Q 

<  P 

q:  -I 

U)  < 


A? 


VELOCITY  360.  720  KM/HR 


Z2  '^HAD/KII^  J  '  GRAVITY  FIELD  M(M)ELS 


O  52 

°  3 


m  >* 

,4  c> 

»  ^ 
Q  w 
o  o 
s  o* 
» 
H  05 
Eh  fc 
O 

O  >  C3 

O  <  o 

.  .  05  ►J 

rH  U 

i 

D  W 


9  w  w 

O  ^  M 
?  Q 
7  m  < 
‘  05 

W  C5 
Eh 

05  N 
O  N 

Q  < 

<  O 

_  M 

o  O  CQ 

®  M  flZ 

O  H  w 


o 

o 

o 

o 

o 

o 

o 

o 

o  >  aQ 

o 

o 

•o 

o 

O  £h 

» 

• 

• 

' 

<  isj 

to 

02 

rH 

O  15 

320  KM 


BASELINE  CONDITIONS 


GYRO  RANDOM  DRIFT 
(DEG/HR  RMS) 


VERTICAL  POSITION 
ERROR  (FT  RIMS) 


lU 


H 

£ 

Q 


kl 


^  I 

I 

o  is 

ut  so 

5  iS£ 

S  oig 

««  oc  j 

<  >•  UJ  I 

cal  o>n 


noise  (E^/RAD/SEC) 


50 


Summary  of  Parametric  Study  Conditions 
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GRADIOMETER 
WHITE  NOISE 
(E.2/RAD/S) 


Navigation  Performance  Analysis 
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TERRAIN  BASED  GRAVITY  MAP 
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TKItRAIN  DERIVED  MAP  PERFORMANCE 
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OCEAN  GRAVITY  IMP 


61 


Navigation  Performance  Anuiysis 

(Position) 
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FOREWORD 

This  docLiment  contains  materiai  used  in  a  presentation 
given  by  The  Analytic  Sciences  Corporation.  The  material  is  not  in¬ 
tended  to  be  self-explanatory,  but  rather  should  be  considered  in  the 
context  of  the  overall  presentation. 
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Summary  findings  from  the  rail  tests 
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DATA  COLLECTION  SCENARIO 
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large  variability  in  elevation 

•  Five  Dodge  tracks  and  three  Albuquerque 
tracks  provided  by  Bell  Aerospace.  Inc. 


TYPICAL  NAVIGATION  DATA  (STAGE  I  OUTPUTS) 

FROM  DODGE  ROUTE 
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Data  are  from  track  6  (Dodge  City  to  La  Junta) 

Position  and  elevation  profiles  reflect  use  of 
transponder  truth  data 

“Corrected”  scales  refer  to  transponder 
smoothing  of  axle  tachometer  output 


TYPICAL  NAVIGATION  DATA  (STAGE  I  OUTPUTS) 

FROM  DODGE  ROUTE  (cont.) 
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“Corrected”  scales  refer  to  transponder 
smoothing  of  axle  tachometer  output 


TYPICAL  NAVIGATION  DATA  (STAGE  I  OUTPUTS) 
FROM  ALBUQUERQUE  ROUTE 
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Data  are  from  track  8  (Albuquerque  to  La  Junta) 

Position  and  elevation  profiles  reflect  use  of 
transponder  truth  data 

“Corrected”  scales  refer  to  transponder 
smoothing  of  axle  tachometer  output 


TYPICAL  NAVIGATION  DATA  (STAGE  I  OUTPUTS) 
FROM  ALBUQUERQUE  ROUTE  (cont.) 


(Bep)  tsBuy  ||0U  uuoptid 


(Sap)  Moud  p3)03JJ03 


74 


Data  are  from  track  8  (Albuquerque  to  La  Junta) 

Pitch  and  roil  magnitude  are  generally  about 
twice  as  large  as  those  for  Dodge  route  ( <  4deg) 

“Corrected”  scales  refer  to  transponder 
smoothing  of  axle  tachometer  output 
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ASSESSMENT  OF  GGSS 
NAVIGATION  SYSTEM  PERFORMANCE 


o 

CL 

'c 

o 

T3 

(fi 

ca 

o 


tL 

O 

•  WM 

o 

o 

o 

> 

•S 


CB 

O 

ID 

D 


O) 

■2-2  <5 

D)  — 

Q  5  •“ 

•o  'O  O 

5  <D  2 
E 

^  S  i 

E  (5  y 
S  w  S 

af  ® 


^3 

|t 

•=  o 

■o  S 

©  E 
to  o 

s 

to 

0  ^ 
O 

O)  o 
E 

V) 


(0 


o  ^ 

r* 

O  55 

.p  O 
CO  a 
>  _ 
0  ^ 

®  s 

«  CO 
0  O 

TD 
D 

’O)  i 
o  -D 

_  C 


D> 

C 

o  "D  -D 

0* 

■D 
D 


0 

•D 

D 


O 

a 
(/> 
c 

M  to 

0  k. 
— 

TD  U> 

B  £ 

o  w 

0  D 

k. 

k-  ~ 

O  0 

o  m 


•  • 


0 

0 


0 

0 

3 

O’ 

c 

o 

0 

c 

0 

E 

0 

0 

0 

0 

0 

< 


8  § 

§  S- 

i  ^ 

*o  > 

V>  'k. 

O  tl) 

to  '9 
tl) 

g>  ^ 

o  E 

0  O 

9  o 
E  to 
0 

0  ^ 

.i  S 

4-* 

8  = 

w  w 

®  T3 
>  ® 

4^^  0 

S  - 

S*  o 

D  O 
O  0 

s  « 
>.  ^ 

^  i 

c  o 

<  o 


76 


Compare  corrected  positions  with  corresponding  transponder  values 
Compare  total  distance  covered  between  first  and  last  common  transponders 

Analyze  position  profiles  along  straight  stretches  (identified  from  topographic 
maps) 
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COMPARISON  OF  CORRECTED  POSITIONS  AND 
APPARENT  TRACK-TO-TRACK 
DIFFERENCES  BETWEEN  TRANSPONDER  LOCATIO 
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outputs,  these  track-to-track  diffei 
temporary  cietours  to  track  sidings 


TYPICAL  COMPARISON  OF  CORRECTED  POSITIONS 

AND  TOPOGRAPHIC  MAPS 

#  18-kin  Straight  section  identified  between  ",  101  ‘‘02'07")  and 

(37“57'25",  lOO^SO  M") 
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#  No  significant  deviations  from  straightness  are  apparent 

#  implies  that  positions  between  transponders  are  sufficiently 
accurate  to  support  registration  of  gravity  quantities 


SUMMARY  OF  NAVIGATION 

ASSESSMENT 
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between  transponders  for  registration  of  gravity 
quantities 


T-20503 
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ANALYSIS  OBJECTIVES 
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Estimate  errors  caused  by  seif-gradients  and 
apply  compensations 


TIME-SERIES  DATA 
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Acceleration  in  piatform  frame,  three  components 
(x,  y,  z)  sampied  at  1 28  Hz 
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DISTRIBUTION  OF  VARIANCES  FOR  TRACK  11, 

DODGE  CITY 
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which  is  predominantiy  in  the  aiong-track  direction 

Anisotropic  variances  are  result  of  correlation  between  rail-track 
and  gravity-field  geometries 


DISTRIBUTION  OF  VARIANCES  FOR  TRACK 

DODGE  CITY 


which  is  predominantiy  in  the  aiong-track  direction 

Anisotropic  variances  are  result  of  correlation  between  rail-track 
and  gravity-field  geometries 


OVERVIEW  OF  ALBUQUERQUE  DATA 
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Consistency  from  track  to  track 


POTENTIAL  BENEFIT  OF  SMOOTH 
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Foregoing  analysis  needs  refinement  but  generai  magnitude  of 
benefit  gained  by  accounting  for  anistropy  is  quite  encouraging 


DISTRIBUTION  OF  VARIANCES  FOR  TRACK  8, 

ALBUQUERQUE 
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variances  (54.5%);  this  is  consistent  with  typicai  gravity-fieid  modeis 

NN  and  EE  variances  are  similar,  because  along-track  direction 
changes  significantly  along  Albuquerque  tracks;  this  differs 
from  Dodge  City  tracks 


POWER  SPECTRA  FOR  STATIONARY 

AND  MOVING  GGSS 


«0  QO 
I 


(zH/za)  asd 


90 


STA-11  and  STA-12  data  yield  self-noise  power  spectra  (PSDs) 
before  and  after  tracks  11  and  12  were  traversed 

Low-noise  process  appears  in  3TA-12 

GGSS  iowpass  filter  roils-off  spectra  above  0.01  Hz 


POWER  SPECTRA  FOR  STATIONARY 
AND  MOVING  GGSS  (Cont.) 
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Increased  spectrum  level  In  TRK 12  at  low  frequencies  supports 
hypothesis  that  seif-noise  increased  during  track  12  traverse 


ENVIRONMENTALLY  INDUCED  NOISE 
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i£3)  )U3!O«J0  A^iABiO 


(3)  tuatpBJQ  A«abjo 


Distance  (km) 


OBJECTIVE  MEASURE  OF  REPEATABILITY 
BETWEEN  TRACKS  8  AND  10 

Squared  Coherences,  Tracks  8  and  10 
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Squared  coherence  is  estimated  using  a  canonical-variates  ^tate-spac 
modeling  technique 

NN,  EE,  and  DD  have  largest  coherences  (C),  because  their  signai-tc 
noise  ratios  (SNR)  are  largest  (  SNR  ) 


EXAMPLES  OF  TIME  SERIES  HAVING 

HIGH  COHERENCES 
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Lowpass  filter  has  half-power  frequency 
coressponding  to  60-km  wavelength 


EXAMPLES  OF  TIME  SERIES  HAVING 

LOW  COHERENCES 
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Lowpass  filter  has  half-power  frequency 
coresponding  to  60-km  wavelength 


REPEATABILITY  BETWEEN  TRACKS  11  AND  12 

Squared  Coherences 
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Gradients  involving  gravity  changes  in  east  direction  (predominantly 
along-track  direction)  have  small  signal  levels  and  are  incoherent 
(except  for  NE  at  wavelengths  longer  than  25  km) 


ACCELERATION-INDUCED  GGSS 
MEASUREMENT  ERROR 

Squared  Coherence  Between  GGi1,  inline,  and  Acceteration  Power,  Track  11 
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Note:  inline  signals  near  7x1 0~^  Hz  are  hetrodyned  to  low  frequencies  when  the  carouselied 
GGI  outputs  are  transformed  to  local-level  coordinates 


EXAMPLE  OF  ACCELERATION  COMPENSATION 
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Example  shows  that  compensation  is  broadband 


COMPENSATION  VARIANCES 

Variances:  GGSS  Data  and  Compensation  Signals,  Track  1 1 
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gradient  elements 

Largest  fractional  improvement  occurs  for  NE,  ND,  EE,  and 
ED  gradients 


COMPENSATION  POWER  SPECTRA 

Spectra  for  ED,  Track  1 1 
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Apparent  over  compensation  at  low  frequencies  may  be  spectrum 
estimation  limitation  with  finite  time  series 

Conclusion:  Additional  vibration  compensation  technique  should 
be  considered 


ENROUTE  SELF-GRADIENT 


102 


•  unginai  motivation:  rail  car  gradiometer  calibration 

•  Seif-gradients  are  biases  in  body  frame 

•  DIrectionai  changes  of  survey  vehicie  moduiate  seif-gradients 
In  local-level  frame 

•  Transform  measurements  into  body  frame  and  subtract  a  mean 


VEHICLE  SELF  GRADIENTS 
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Thus  we  are  motivated  to  transform  the  measurements  into  a  frame  in 
which  the  gradiometer  platform  does  not  rotate 


VEHICLE  SELF  GRADIENTS  IN 
THE  BODY  FRAME 
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Foregoing  formulation  considers  only  second-order  gradients  which 
transform  as  above  —  higher-order  self-gradients  have  smaller  effect 


TYPICAL  ALBUQUERQUE  TRACK 
SELF-GRADIENT  COMPENSATION 


100  150  200  250  300  350  400  450  500 

Aiong-Track  Distance,  Track  8  (km) 


TYPICAL  ALBUQUERQUE  TRACK 
SELF-GRADIENT  COMPENSATIQN 
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Normalized  by  the  standard  deviation  of  the  corresponding  gradient  on  the  outbound  track 


IMPLICATIONS  OF  SELF-GRADIENT 
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Similar,  but  smaller  effect  on  the  Dodge  City  tracks 


PRIMARY  CONCLUSIONS  FROM 

DATA  ANALYSIS 
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but  they  account  for  onfy  a  fraction  of  the  total  measurement  noise 
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RATIONALE  FOR  OBSERVED  DIRECTIONAL 
BEHAVIOR  OF  GRADIENTS 
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rms  M.-  gradient  is  large  due  to  mass  contrast  on  opposite  sides  of  railroad  grade 

rms  4  X  -  gradient  Is  moderate  for  same  reasons  as  above  but  offset  by  reduced  rms  of  cross 
gradients  vs  Inline  gradients 

rms  ^lZ‘  gradient  Is  small  due  to  cross-track  uniformity  of  roadbed 
rms  z  z  -  gradient  behaves  like  ^-gradient  (by  Laplace’s  Equation) 


SUMMARY  OF  DIRECTIONAL  DEPENDENCE  OF 
HORIZONTAL-INLINE  GRADIENTS 
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Also  note  higher  rms  on  Albuquerque  route 

ou(n,  12) 


IDEALIZED  MODEL  OF  GRADIENT  ANISOTROPY 
ALONG-TRACK  VS  CROSS-TRACK 
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QUICK  LOOK  AT  EFFECT  OF  GRADIENT  ANISOTROPY 
ON  RAILBED  DEFLECTION  SURVEY  DENSIFICATION 


o  Ov 
n  30 

•-  I 


o 

••• 

M 

C 

(0 


o 

Ui 


c 

o 

(Q 

O 

O 


0)  o 

E  o 

0)  O) 
•C  c 
o 

C  _ 

o  O 

(3  c 

S  o 

«  ? 

C  o 

a  r* 

X)  0. 


c 

o 


(0 


o 
n 

M  d 
in 

o 

2"  ^ 

w  .2 

0)  ^ 
a  8  ^ 

1  'S 

</)  T3 


■g 

o 

q: 


O  0. 


0) 

c 

o 

•«M 

ra 

E 

in 

<u 

c 

o 

•«w 

o 

o 

o 

•a 

in 

E 


C? 

o 

o 

c 

<0 

in 

c 

o 


<0 

0)  < 


¥0  € 


ka 

o 

cST' 

o 

o 

O 

x: 

II 

♦i* 

o 

w 

XJ 

a 

2 

k» 

vMji 

CB 

CL 

E 

o 

o 

7b 

E 

£i 

CO 

C 

< 

e 

> 

0) 

£ 

CO 

!««• 

> 

2 

O) 

« 

c 

o 


T3 

■o 

(Q 

•o 

c 

ra 

o 

13 

o 

E 


> 

2 

Oi 

75 

o 

35 

2 

CO 

(0 


CB  g 

to  I 

2  §• 

0)  »- 

•c  » 

c  I 

§.  5 

3  tS 
«  C 

c  2 

O  (M 

o  Q 
*0 

•o 

1  o^S 

2  o 

«  « 
f~  C 
>  0) 
>  (0 


■D 

C 

(B 

Q 


II  S 
\  in  ^ 

a  c  VI 

■O  S  «J 

c  S  S 

c  ®  ^ 

3  C 

in  *0 

CO  S 

(B  Q. 

■D  ® 

®  o 

t 

2  ® 

1 

2  w 


0> 

(0 

(B 

O 


w 

«> 

o 

c 

(0 

in 

'o 

I. 

2 

CO 

c 

CB 

%m 

E 

3 

E 

1 

E 

CO 

2 

vm» 

a 

E 


1 

a 

JUJ> 

rH 

L 

Cl 

VI 

< 


o 

o 

Q. 

35 

3 

E 

(0 

o= 

•k 

JIB 

u 

CB 

k. 

I 

O) 

c 

o 

(B 

B> 

8  >5 

2 

in  w 


c 

O  ^ 
%  ® 


o 

o 

c 

o 

3S 

o 

2 

•3 

O 

o 


«> 

CO 

IB 

o 


113 


Implies  maximum  transfer  distance  increases  to 

A'  <  D;i  r'l^max/o-^  =  A /a) 
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ASTRO  SITE  LOCATIONS 

TRUTH  DATA 
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COMPARISONS  WITH  TRUTH  DATA 
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Regional  Mean/Trend  Removed 


TRUTH  DATA  VS  GGSS  ESTIMATES 
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(Albuquerque,  NM)  (La  Junta,  CO)  (Dodge  City,  KS) 

Astro  Station  Distance  along-track  (full  scale  =  880  km) 


SUMMARY  FINDINGS  FROM 

RAIL  TESTS 
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Aiong-track  smoothness  will  reduce  railroad  astro  survey 
densification  costs  significantly 


Obtaining  Harth  Surface  and  Spatial  Deflections  of  the  Vertical  frQm 
Free-Air  Gravity  Anomaly  and  Elevation. Data  WithQjJt  Density J\ssumPtiQn.s 


DAVID  M.  GLEASON 

GEOPHYSICS  LABORATORY 
HANSCOM  AFB,  BEDFORD,  MA.  01731 


ABSTRACT:  Moritz  (1980)  presents  a  density-free  scheme  allowing  for  the 
analytical  continuation  of  a  given  set  of  free-air  gravity  anomalies  to  any 
desired  level  surface  if  a  corresponding  set  of  elevations  (e.g.,  above  MSL) 
is  available.  An  efficient  spectral  implementation  of  this  scheme  is 
discussed  by  Sideris  (1987).  A  subsequent  spectral  execution  of  the 
planar  Vening-Meinez  equation  on  the  continued  anomalies  yields 
deflections  of  the  vertical  on  the  chosen  level  surface.  The  deflections 
are  brought  back  to  the  Earth's  surface  via  a  spectrally  implemented 
Taylor  series.  The  series'  convergence  rate  depends  on  a)the  ruggedness  of 
the  local  topography  and  b)the  resolution  of  the  input  gravity  and  elevation 
gridded  data.  Deflections  at  a  constant  altitude  above  the  level  surface 
are  obtained  through  a  routine  spectral  execution  of  the  planar  upward 
continuation  integral.  Two  sites,  having  diverse  topographies,  were 
surveyed  for  1'  by  1'  mean  free-air  anomaly  and  elevation  values  and  for 
smaller  sets  of  astronomically-determined  deflections  to  serve  as  control 
or  "truth"  values.  In  a  topographically-tranquil,  but  gravimetrically 
turbulent  Oklahoma  site  the  overall  RMS  of  the  differences  between  true 
and  predicted  deflections  was  0.3  arc  secs  and  in  a  rugged  New  Mexico 
site,  using  less  reliable  truth  data,  the  RMS  was  0.6  arc  secs.  Potential 
pitfalls  of  the  2  dimensional  Fast  Fourier  Transform  pair  are  discussed 
with  an  emphasis  on  unwanted  circular  convolution  effects  which,  if 
unaccounted  for,  can  increase  the  error  in  predicted  deflections  by  as 
much  as  100%. 
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OBTAINING  EARTH  SURFACE  AND 
SPATIAL  DEFLECTIONS  OF  THE 
VERTICAL  FROM  MEAN  FREE-AIR 
GRAVITY  ANOMALY  AND  ELEVATION 
DATA  WITHOUT  DENSITY 
ASSUMPTIONS. 

DAVID  M.  GLEASON 
GEOPHYSICS  LABORATORY 
HANSCOM  AFB,  MA.  01 731 
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p 


Range's  Thm.  states  one  can  always  find  a 
harmonic  function  T  ,  arbitrarily  close  to 
"^EXTERNAL’  regularly  continued 

(be  it  upward  or  downward)  from  the  ground 
to  a  chosen  level  surface. 
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So 


I 

3Ag'p, 

^9p  =  p  +  ^p'  02 


+ 


•  The  Ag'  set  reflects  the  earth's  exterior 

(not,  interior)  gravity  field.  Used  in  StoKes' 
formula,  it  yields  a  T'  which  is  harmonic  above 
the  chosen  surface  and  which  agrees  with  the 
actual  T  on  and  above  the  Earth's  surface. 


•  So,  under  such  a  continuation  of  T^xj, 

masses  outside  the  chosen  surface  are,  in 
.effect,  shifted  to  its  interior. 
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Again 


Ag'pt  =  g’^p+  g^p'  +  g^p-  +  g^p'  + 
where 

g^p  =  the  observed  Agp 


2 


g 


p’ 


Etc.,  Etc. 
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•  In  g''p-  =  -zp-(9Ag/3z),  the  partial  will 
be  treated  as  a  planar  surface  operator. 


i.e. 


I 

=  -7  • - 

P  dz 


=  -z 


r.  g  (x.y)  -  g  (Xp.Vp) 

pJJ - 2 - 23^dxdy 

■“  [(x-Xp)  +(y-yp)  ] 


which  is  a  convolution. 
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NOTES: 


1) .  Due  to  only  a  finite  grid  of  Ag  input  values, 
the  spectrum  of 

^gREDUCED  =  ^gOIVEN  '  ^gS.H.  EXPAND 
is  computed . 

2) .  As  in  all  applications  of  the  2D  EFT  pair, 
be  aware  of 

•  Aliasing 

•  Spectral  Leakage 

•  Circular  (non-linear)  Convolution  Effects 
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NOTES:  , 

1.  After  obtaining  the  gridded  set  of  Ag' 
values  on  the  chosen  surface,  one  can 

immediately  obtain  gridded  sets  of  and  Tj' 
deflections  on  the  chosen  surface  via  a 
routine  spectral  execution  of  the  planar 
Vening  Meinez  eqn  using  the  applicable 
transfer  functions. 

2.  One  can  then  obtain  gridded  sets  of 
spatial  deflections  at  a  constant  altitude  h 
above  the  chosen  level  surface  through  a 
routine  spectral  execution  of  the  planar 
upward  continuation  integral  (using  the  u.c. 

transfer  function  e'®*^). 

3.  One  can  efficiently  obtain  gridded 
deflections  on  the  irregular  Earth  surface 
via  a  spectral  execution  of  the  Taylor 

Series  linking  ^'p.  and  ii'p.  to  ^p  and  tip  . 
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Topography  of  Central  3°  by  3°  New  Mexico 
Area. 


gn-1  Values 


g  MAX  = 
86  mgals 


*  f  'V^'’ 


gn-2  Values 

g^MAX  = 

72  mgals 
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gn=2  Values 


g^MAX 
72  mgals 


gn-3  Values 

g^MAX  = 

60  mgals 


Overall  RMS  values  of  378  (Astro-Predicted) 
^  and  Ti  differences,  in  arc  secs,  from  1 '  New 
Mexico  data  and  the  Newer  Astro  set. 


n  Truncation  Level 

Tl 

0 

CD 

• 

O 

0.75" 

1 

0.59" 

0.74" 

2 

0.59" 

0.74" 

3 

0.59" 

0.74" 

NOTE: 

^Predicted  =  ^  while 

^Astro  “  <I>  -  (j) 

where  (j)  -  (])  =  (f  h/R)sin2(j),  is  the 

well-known  reduction  for  the  normal  curvature 
of  the  plumb  line. 


•  The  predicted  (gridded)  deflections  were 
interpolated  to  the  astro  locations. 
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From  Data  Types  and  Their  Spectral 
Properties,  by  K.P.  Schwarz: 


Percentage  of  Total  Value,  by  Harmonic 
Degree  n,  for  T,  and  T,,  (i.e.,  Ag  and 

3  Ag/3z) . 


low 

medium 

high 

very  high 

n  e  (2,36) 

(37,360) 

(361 ,3600) 

(3601 ,36000) 

(5°  grid) 

(30'grid) 

(3-  grid) 

(18"  grid) 

Tz  22.5 

41.9 

32.7 

2.8 

o 

• 

o 

k-T 

0.8 

39.0 

60.2 

>  NOTES: 


A 1' grid  <=>  nmax=  10800 
A  30"  grid  <=>  nmax=21600 
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SUMMARY: 


1 .  The  spectral  approach  allows  for  efficient 
predictions  of  deflections  and  height 
anomalies  at  a  resolution  matching  the  input 
data. 

2.  The  n=1  topographic  corrections  were 
beneficial  in  near-mountain  splotches  where 

the  input  Ag  and  h  data  was  reliable  but  were 
detrimental  in  such  areas  where  the  data  was 
suspect. 

3.  The  extraction  of  ultra-high  frequency 
information  from  lower  order  1'  or  30"  mean 
gradients  is  questionable. 

4.  Noise  in  such  input  data  might  render 

higher  order  (e.g.  g3,4,5,6 . j  corrections 

meaningless. 

5.  Interpolated  grids  of  anomalies  and 
heights  fail  to  account  for  higher  frequency 
terrain  effects. 
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Distinguishing  Nuclear-  from  Conventionally-  Armed 
Cruise  Missiles  with  a  Gravity  Gradiometer 


Dr.  John  A.  Parmentoia 
Center  for  Science  and  International  Affairs 
John  F.  Kennedy  School  of  Government 
79  John  F.  Kennedy  Street 
Harvard  University 
Cambridge,  MA  02129 


Abstract 


I  will  discuss  an  analysis  of  an  application  of  the  gravity 
gradiometer,  that  has  been  designed  at  Draper  Laboratories,  which 
might  be  useful  during  missile  production  for  distinguishing 
conventional  from  nuclear  armed  cruise  missiles  in  a  nonintrusive 
way.  The  motivations  for  exploring  this  potentially  important 
application  of  this  device  will  also  be  discussed. 
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Distinguishing  Nuclear  From 
Conventionally  Armed  Cruise  Missiles 

I.  Brief  Introduction 

Arms  control  and  START  : 

Long-range  nuclear  sea-launched 
cruise  missiles  or  SLCMs. 

For  the  U.S.  these  are  the  various 
versions  of  the  Tomahawk. 

For  the  Soviets  these  are  the  SS-N- 
21s. 

Cruise  Missiles  : 

Pilotless  aircraft 
Sophisticated  autonomous  guidance 
Fly  low  to  avoid  radar 
Nuclear  or  conventional  payloads 
Variety  of  launch  platforms 
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Cruise  missiles  along  with  their 
protective  cannisters  are  loaded 
into  launchers. 

Upon  firing  the  cruise  missile 
breaks  through  the  cannister  while 
being  propelled  by  a  rocket  motor. 

Rocket  motor  burns  until  cruising 
speed  is  reached,  then  a  turbofan 
engine  takes  over. 

Guided  to  its  target  by  local  terrain 
maps  and  for  the  conventional 
version  an  additional  digital  image 
of  target. 

The  issues  I  want  to  address  in  this 
talk  are  : 

(1)  Why  are  long-range  nuclear  SLCMs 
an  important  issue  ? 

(2)  What  problems  do  SLCMs  pose  for 
verification  ? 

(3)  How  could  a  gravity  gradiometer 
contribute  to  a  SLCM  verification  scheme  ? 
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The  quick  answers  to  these  questions 
are  : 

(1)  Because  long-range  nuclear  SLCMs 
are  a  contentious  issue  within  START. 
There  are  important  political,  military,  and 
economic  factors  which  contribute  to  this 
situation. 

(2)  Several  design  characteristics  of 
long-range  SLCMs  make  verification  of 
limits  on  SLCMs  relatively  hard,  but  far 
from  impossible  in  my  opinion. 

(3)  Because  long-range  nuclear-  and 
conventionally-armed  SLCMs  have 
significantly  different  internal  mass 
distributions  a  gravity  gradiometer  could 
effectively  distinguish  between  them.  This 
would  ensure  that  nuclear  cruise  missiles 
are  not  being  falsely  counted  as 
conventional  under  an  arms  control  treaty. 
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II.  Why  worry  about  Sea-Launched  Cruise 
Missiles  or  SLCMs  ? 

There  are  military,  economic,  and 
political  dimensions  to  this  question. 

Military: 

Sneak  attack  of  bomber  bases  and 
command  and  controi  system  with  a 
small  number  of  stealthy  SLCMs.  No 
early  warning  ! 

Economic: 

The  costs  of  deploying  an  effective 
early  warning  system  are  high. 

Political: 

SLCMs  contributed  to  the  stopping 
of  START.  The  most  contentious 
aspect  has  been  verification. 

Simplest  situation  is  a  ban  on  both 
conventional  and  nuclear.  U.S.  Navy 
opposes  a  total  ban,  because  it  likes 
the  conventional  anti-ship  version. 
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Domestically  the  public  has 
demanded  reductions  and  it  is 

unlikely  that  SLCMs  will  be 

excluded. 

Arms  control  might  provide  the 

means  for  controlling  the  threat, 
save  the  country  some  money  for  an 
early  warning  system,  and  solve 

some  political  problems. 


138 


III.  What  are  long-range  SLCIUIs  ? 

I  will  restrict  my  discussion  to  U.S. 
cruise  missiles. 

U.S.  Tomahawk: 

Land-attack  nuclear  has  a  range  of 
2800  km. 

Land-attack  conventional  has  a 
range  of  1500  km.  Has  shorter 
range  due  to  much  longer  and 
heavier  warhead,  less  fuel,  and 
more  guidance  equipment  for 
targeting. 

Anti-ship  conventional  has  a  range 
of  500  km.  Even  more  guidance 
equipment,  because  it  goes  after 
moving  targets  and  has  less 
efficient  engines. 

Verification  probiems: 

Nuclear  and  conventional  versions 
aii  have  the  same  airframe. 
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Minor  visual  external  differences 
are  not  useful  for  distinguishing 
between  them. 

Launchers  are  dual-purpose. 

This  impiies  that  NTM  is  useless.  More 
intrusive  forms  of  monitoring  are  required 
such  as  on-site  inspections  of  : 

Storage  facilities 

Service  facilities 

Testing 

Depioyment  modes 

Production  -  Gravity  Gradiometer 
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IV.  Some  difficulties  with  SLCM 
verification 

Simpiest  situation: 

Total  ban  on  all  long-range  SLCMs  - 
Requires  the  dismantling  of 
infrastructure  needed  to  produce, 
store,  service,  test,  and  train  with 
these  weapons. 


Much  harder  situation  : 

Problems  occur  when  a  category  is 
allowed.  Infrastructure  for 
conventional  anti-ship  version  of 
Tomahawk  usefui  for  producing 
long-range  nuclear  SLCMs. 

Difficulties  with  verifying  a  ban  on 
nuciear  with  a  iimit  on  conventionai  : 
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Possibility  that  nuclear  warhead  is 
mated  to  a  cruise  missile  in  the 
factory  and  designated  as 
conventional.  This  could  be  detected 
nonintrusively  with  nuclear 

radiation  detectors  or  intrusively 
with  x-rays  or  a  beam  of  neutrons. 

It  is  not  necessary  that  a  nuclear 
warhead  be  mated  in  the  factory  to 
cheat.  What  is  needed  is  a  tested 
production  line.  Place  a  dummy 

warhead  with  the  same  mass 
distribution  as  nuclear  warhead. 
This  could  be  detected 

nonintrusively  with  a  gravity 
gradiometer. 


If  an  effective  nonintrusive  method  of 
distinguishing  nuclear  and  conventional  can 
be  found  then  tags  can  be  used  to  identify 
legal  cruise  missiles  throughout  their  life- 
cycle. 
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V.  Using  the  Gravity  Gradiometer  for 
Production  Monitoring 

Mass  density  distributions  of  a  nuclear 
and  conventional  Tomahawk  : 

Nuclear  warhead  is  more  than 
twice  as  dense,  half  as  long,  and 
more  forward  in  location  than  the 
conventional  warhead. 

More  than  half  the  volume  of  fuel  in 
the  nuclear  version  is  in  roughly  the 
same  location  as  the  warhead  in  the 
conventional  version;  however  the 
density  of  the  fuel  is  roughly  20  % 
less  than  the  warhead. 


Draper  three  gradiometer  device: 

A  properly  packaged  and  field 
tested  device  would  cost  about  2 
million  dollars,  i.e.  the  price  of  a 
Tomahawk. 
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This  version  wouid  minimize  errors 
due  to  jitter. 

Accurate  to  the  1  Eotvos  unit  level. 
Response  time  is  relatively  short. 

Place  on  a  tripod  in  a  room. 

A  Tomahawk  would  pass  by  the 
device  at  a  certain  rate  while  the 
gradiometer  responded  and  provided 
a  complete  scan  of  the  missile. 

It  would  be  desirable  to  allow  for 
several  different  scans  along  the 
missile  length  at  different  distances 
from  the  missile  axis. 

Since  the  response  time  of  the 
device  is  relatively  short  the  time 
required  to  scan  a  cruise  missile  is 
not  a  significant  factor  in  the 
monitoring  process. 


The  resolution  of  the  device,  i.e.  its 
capability  of  discerning  details 
about  the  individual  components  of  a 
cruise  missile,  is  easy  to  control  by 
limiting  how  close  to  the  missile  the 
gradiometer  can  approach. 

Estimating  the  distance  for  a 
hypothetical  measurement: 

Specify  the  desired  level  of 
accuracy.  Assume  for  the  moment 
that  cruise  missile  is  a  very  long 
cylinder  of  uniform  density. 

An  inaccuracy  in  a  measurement  of 
the  radial  gradient  of  the 
gravitationai  field  translates  Into 
an  uncertainty  in  the  determination 
of  the  mass  density  or 
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8F'/F  =  8p/p 

dp/p  =  r2  dF/[  2TTa2  p  G  ] 

<5F'  =  1  Eotvos  unit 

dp/p  =  0.02 

a  =  0.265m 
p  =  2  X  10^  kg/m^ 
r  «  1m 

Model  of  mass  density  used  in  the 
caicuiations  : 
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The  aluminum  skin  of  the 
hemispherical  nose  and  Its  volume 
mass  density  are  approximated  by 
points  located  at  their  respective 
centers  of  mass. 

The  remaining  components  including 
the  airframe  are  treated  as  lines  of 
mass  along  the  axis  of  the  missile. 

For  this  model  the  gradiometer  will 
be  sensitive  to  the  mass  per  unit 
length. 

Results  : 

Simulation  of  the  radial  gradient  of 
the  radial  component  of  the 
gravitational  field,  Frr,  produced 
by  a  conventional  and  nuclear  cruise 
missile  along  their  length. 

If  measurements  are  not  made 
continuously  along  the  missile  axis, 
then  the  stepsize,  i.e.  the  distance 
between  measurements,  must  be 
less  than  the  size  of  an  object  in 
order  to  resolve  it. 
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eotvos  units 


Radial  Gradient  Comparison  at  r  =  0.5m 
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•otvo*  unit* 


Radial  Gradient  Comparison  at  r  3  1m 


Z  (  m  ) 
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•otvos  units 


Radial  Gradient  Comparison  at  r  =  0.5m 


Z  {  m  ) 
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Model  of  Conventional  SLCM 


Section 

Component 

Mass 

Length 

Radius 

Skin  Thickness 

Average  Density 

1  * 

Guidance  System 

68  kg. 

.646m 

.252m 

.013m 

.61gm/cc 

2 

Fuel 

27 

.128 

.252 

.013 

1.07 

3 

Warhead 

456 

1.770 

.252 

.013 

1.30 

4** 

Fuel 

176 

1.400 

.244 

.021 

.67 

5 

Engine 

59 

1.640 

.252 

.013 

.18 

6a 

Fuel 

178 

.652 

.220 

.000 

1.80 

6b 

Rocket 

121 

.652 

.265 

.000 

2.70 

7 

Skin  or  Airframe 

365 

5.590 

2.70 

r  •  Nose  is  assumed  to  be  a  hemisphere  of  radius  .252m. 

4**  -  Wings  are  inciuded  in  airframe  which  accounts  for  smaiier  inner  radius.  Top 
surface  of  wings  is  assumed  to  have  an  area  of  1 .02  m^  and  the  mass  of  the  wings  is 
assumed  to  be  52.5  kg. 
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Model  of  Nuclear  SLCM 


Section 

Component 

Mass 

Length 

Radius 

Skin  Thickness 

Average  Density 

1  * 

Guidance  System 

46  kg. 

.458m 

.252m 

.013m 

.61  gm/cc 

2a 

Warhead 

123 

.866 

.130 

.000 

2.70 

2b 

Fuel 

123 

.866 

.252 

.013 

.97 

3 

Fuel 

260 

1.220 

.252 

.013 

1.07 

4  *  * 

Fuel 

176 

1.400 

.244 

.021 

.67 

5 

Engine 

59 

1.640 

.252 

.013 

.18 

6a 

Fuel 

178 

.652 

.220 

.000 

1.80 

6b 

Rocket 

121 

.652 

.265 

.000 

2.70 

7 

Skin  or  Airframe 

365 

5.590 

2.70 

r  -  Nose  is  assumed  to  be  a  hemisphere  of  radius  .252m. 

4**  •  Wings  are  i.ncluded  in  airframe  which  accounts  for  smaller  inner  radius.  Top 

surface  of  wings  is  assumed  to  have  an  area  of  1 .02  m^  and  the  mass  of  the  wings  is 
assumed  to  be  52.5  kg. 
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17th  Gravity  Gradiometry  Conference 

12  -  13  October,  1989 

Air  Fo'”ce  Geophysics  Laboratory 
Hanscom  AFB 


Advances  in  Dynamic  Estimation 


Dave  Sonnabend 
Jet  Propulsion  Laboratory 
California  Institute  of  Technology 

Abstract 


In  prior  years  I  have  talked  about  magnetic  isolation  of 
instruments,  with  only  short  allusions  to  our  work  in 
dynamic  estimation  to  deal  with  rotation  correction  in 
floated  gradiometers.  This  year's  talk  will  be  almost 
entirely  devoted  to  estimation.  As  the  theory  has  been 
exposed  at  other  conferences  and  seminars,  and  is  a  central 
topic  in  my  book  on  gradiometry,  it  will  only  be  sketched 
here.  However,  there  are  several  new  developments, 
including  improvements  to  our  models  and  filters, 
application  to  a  Lunar  Observer  mission,  and  computational 
techniques  for  dealing  with  self  gravity.  Also,  if  no  one 
from  NASA  Hq.  shows  up.  I'll  discuss  NASA's  latest  plans  in 
gravity  measurements. 
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DYNAMIC  ESTIMATION 
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Measurement  Equations 


ESTIMATION  STATE 
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GRADIENT  SIGNAL  STATISTICS 
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PROPELLANT 
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ACCELEROMETER 

ENSEMBLES 


CONFIGURATION 

(0.5  m  EDGE) 

PROOF 

MASS 

TRIANGLE 

F 

SQUARE 

FC 

TETRAHEDRON 

FC 

OCTAHEDRON 

F 

CUBE  FACES 

F 

CUBE  CORNERS 

FC 

NOISE  LEVEL 

SENSITIVE  1.3  X  10-*^  m/sec^  - 

INSENSITIVE  2.5  x  10-1“  m/sec^  - 
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BUT  THE  PROPELLANT  STICKS  TO  THE  WALL 
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IT  WORKS! 
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RESULTS  ON  THE  ESTIMATION  OF  GEOPOTENTIAL  COEFnCBENTS  FROM 
A  SIMULATION  OF  A  SATELLITE  GRAVITY  GRADIOMETER  MISSION 


Srinivas  V.  Bettadpur,  Bob  E.  Schutz,  John  B.  Lmdberg 
Center  for  Space  Research,  University  of  Texas  at  Austin,  Austin,  Tx  78712 


The  NASA  Satellite  Gravity  Gradiomcter  Mission,  dcsi^^ned  to  measure  the  tensor 
of  gradients  of  accelerations  due  to  gravity,  promises  a  substantial  increase  in  the 
knowledge  of  the  fine  scale  features  of  the  gravity  field  of  the  earth.  One  possible  mis¬ 
sion  scenario  consists  of  the  gradiometer  mounted  in  a  satellite  traveling  in  a  polar, 
frozen  perigee,  drag  free  orbit  and  measuring  the  six  components  of  the  tensor  of  gravity 
gradients  in  a  suitable  reference  frame. 

Some  results  are  reported  from  an  initial  simulation  of  die  estimation  of  the  geopo¬ 
tential  coefficients  from  measurements  made  on  such  a  satellite  gradiometer  mission. 
Using  a  small  reference  gravity  field  (18  by  18  subset  of  a  GEMTl  error  model),  the  gra¬ 
diometer  observations  along  a  true  orbit  were  simulated  in  a  geocentric  equatorial  coordi¬ 
nate  frame.  Zero  mean  Gaussian  random  noise  with  different  standard  deviations  were 
added  to  the  simulated  observations.  During  the  estimation  process,  the  observations 
were  modeled  along  a  nominal  orbit  using  Pines*  fully  nonnalized,  nonsingular  formula¬ 
tion.  To  simulate  a  range  of  orbit  accuracies,  the  dinerences  between  the  nominal  and 
true  orbits  were  varied  from  16,  69  and  19  meters  to  16,  69  and  19  cm  in  the  radial, 
transverse  and  normal  directions,  respectively.  The  geopotential  coefficients  were 
estimated  from  a  least  squares  fit  of  the  simulated  gradiometer  data  in  the  presence  of  dif¬ 
ferent  levels  of  observation  noise  and  orbit  errors.  The  estimated  coefficients  were  then 
compared  to  the  coefficients  of  the  reference  gravity  field,  in  the  sense  of  degree  aver¬ 
aged  errors  and  the  errors  produced  in  a  global  geoid. 

The  results  obtained  from  the  initial  simulations  indicate  that  to  recover  the  global 
geoid  to  about  a  centimeter  root  mean  square  error,  the  instrument  must  have  a  sensitivity 
of  10*^  E.U.  and  the  radial  orbit  accuracy  must  be  within  20  cm.  For  example,  with  errors 
of  16,  69  and  19  cm  in  the  radial,  transverse  and  normal  directions,  and  with  KT*  E.U. 
noise,  the  global  geoid  error  was  0.6  cm  (RMS).  On  the  other  hand,  with  the  same  orbit 
error,  but  with  l(P  E.U.  noise,  the  global  geoid  error  increased  to  45  cm  (RMS). 

Significant  errors  in  the  estimated  coefficients  are  seen  to  be  caused  by  the  adjust¬ 
ments  required  to  model  the  systematic  gradient  residual  due  to  the  point  mass  term  p/r. 
The  permissible  radial  orbit  error  is  governed  by  the  ratio  of  this  systematic  residual  gra¬ 
dient  to  the  noise  level. 

These  results,  while  demonstrating  the  role  of  some  error  sources  in  the  process  of 
estimation,  provide  a  baseline  against  which  the  results  of  approximate  methc^s  can  be 
compared. 
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RESULTS  FROM  THE  ESTIMATION  OF  GEOPOTENTIAL  COEFFICIENTS 


FROM  A  SIMULATION  OF  A 
SATELLITE  GRAVITY  GRADIOMETER  MISSION 


Srinivas  Bettadpur 
Bob  E.  Schutz 
John  B.  Lundberg 


Oct.  12, 1989 
Center  for  Space  Research 
The  University  of  Texas  at  Austin 
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SATELLITE  GRAVITY  GRADIOMETER  MISSION 


Measurement  of  spatial  variation  of  acceleration  due 
to  gravity 

Global,  High  resolution  study 


Goals  : 

Determine  high  degree  and  order  ( ~  180 )  geopo¬ 
tential  field 

Applications  : 

*  Precision  Orbit  Determination 

*  Navigation 

*  Oceanography 
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ANALYSIS  OF  DATA 


Measurements  :  Gradients  of  gravity  in  an  instru¬ 
ment  frame 


Data  : 

*  Orientation  of  the  instrument  frame 

*  Angular  velocity  of  the  instrument  frame 

*  Orbit  of  the  satellite  carrying  the  gradiometer 


Unknown  : 

Coefficients  of  the  spherical  harmonic  expansion 
of  the  geopotential 
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ASSUMPTIONS  FOR  THE  SIMULATIONS 


*  4  second  data  sampling  with  10  or  10  ^  E.U. 
noise 

*  Signal  consists  only  of  the  static  geopotential 

*  Gradients  are  available  in  Geocentric,  Earth  fixed 
frame 

*  Orbit  of  the  satellite  is  separately  available 

*  Error  Sources  : 

Orbit  eiTors 

Observation  noise 
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OBSERVATION  MODEL 


S(tk)  =  V[VU(-f(tk)] 

=  Z  [  “nmij  tk  )  )  Cmn 

n,m 

+  Pnmij  (■i"(  tk  )  )  S„„,  ] 


Yk  )  -  Yk  )  +  ^  [  Yk  ^k 

=  (l|j )  X  +  67^  +  ;  k  =  1 , . 
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THE  ESTIMATOR 


y  =  H  (Tn  )  X  +  e 

E[e]  =  0 

ELee”^  ]  =  I 
X  =  (  h’^  H )"' H’’'’ y 
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DESCRIPTION  OF  SIMULATIONS 


Orbits  : 

*  TRUE  :  32  day  ground  track  repeat  period, 
frozen  perigee,  circular,  polar  orbit 

*  BASIC  NOMINAL  :  Fits  true  orbit  with  errors 
(worst  case) 

16, 69, 19  m  Radial,  Transverse,  Normal 

*  OTHER  NOMINAL  :  Obtained  geornetricaUy 
from  TRUE  and  BASIC  NOMINAL  orbits 

:=  *1^  -  T|i  »  0  <  Tii  <  1 

Best  case  fit :  16, 69, 19  cms 
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DESCRIPTION  OF  SIMULATIONS  (contd.) 


True  Field 

*  1 8  by  18  GEMTl  error  model 

Simulated  Observations  : 

*  Along  TRUE  orbit, 

*  from  TRUE  field, 

*  at  4  sec.  intervals,  for  5  days 


Noise  : 

*  additive  N  ( 0 ,  )  noise 

*  for  a=  10"^  and  a=  10"^  E.U. 
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DESCRIPTION  OF  SIMULATION  (contd.) 


Estimated  Field  : 

y  =  V[VU(fT)]  -  V[V(-!i-)] 

*  Compute  partials  on  the  NOMINAL  orbit 

*  Estimate  same  coefficients  as  in  the  TRUE  field 

Normal  Equations  : 

*  Square  Root  Free  Givens’  Rotations 

*  CRAY  X-MP/24  at  UTCHPC 
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DESCRIPTION  OF  SIMULATIONS  (contd.) 


Description  of  errors  : 
*  Fractional  error 


*  Root  mean  square  global  "geoid"  errors 
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THE  POINT  MASS  TERM  ^ 

r 


Maximum  change  of  gradient  for  PMT 

Radial  distance 

(in  meters) 

100 

50 

5 

1 

0.3 

[  A  Gjj  Imax 
(X10-^E.U.) 

1325 

663 

66 

13 

4 

Maximum  change  of  gradient  for  perturbation  field 

360  by  360  field  OSU86F  ( x  lO"^  E.U. ) 

Displacement 

(in  meters) 

Radial 

Transverse 

Normal 

100 

5 

2 

3 

50 

3 

1 

2 

5 

1 

0.3 

0.1 

0.2 

0.05 

0.02 

0.03 

0.3 

0.016 

0.006 

0.01 
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DEGREE  AVERAGED  ERROR 


DEGREE  AVERAGED  FRACTIONAL  ERROR  AT  END  OF  5  DAYS 
GRAVITY  FIELD  18  BY  18  SUBSET  OF  GEMT1 
EFFECTS  OFORBIT  ERRORS  AND  OBSERVATION  NOISE 


CONTOUR  INTERVRL  -  0.200  CM 

MINIMUM  «  -1. 86  MfiXIMUM  «  1. 46 

MEfiN  =  0. 16  RMS  «=  0. 60 


Gam)  DIFFERENCE  BETUEEN  TRUE  RM)  ESTlMmED  FIELD 
ORBH  ERROR  IR.T.N)  •  (16.6a  18)  CM.  l.E-2  E.a  NOISE 


I  .  I  .  I  . 


CONTOUR  INTERVRL  »  0. 100  M 

MINIMUM  =  -1.09  MflXIMUM  » 

MEAN  =  -0. 17  RMS 
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DEGREE  AVERAGED  FRACTIONAL  ERROR  AT  ENH  OF  5  DAYS 
GRAVITY  FIELD  18  BY  18  SUBSET  OF  GEMT1 
EFFECTS  OFORBIT  ERRORS  AND  OBSERVATION  NOISE 


RESIDUAL  GRADIENTS  FROM  ERROR  IN  PMT 


Radial  orbit  error  =16  meters 

Component 

Gxx 

Gxy 

Gxz 

Gyy 

Gyz 

Gzz 

Avg  resid 
(IQ-^E-U.) 

187 

114 

196 

156 

170 

273 

*  Systematic  residual  gradient  is  compensated  by 
estimated  coefficients 


*  Low  degree  and  order  coefficients  absorb  the  resi¬ 
dual 
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CONCLUSIONS 


*  Residual  gradient  due  to  PMT  affects  errors  in  all 
coefficients. 

*  Allowable  radial  orbit  error  determined  by  ratio  of 
residual  gradient  to  the  noise  level. 

Explicitly  model  the  PMT  residual 
Radial  position  from  GPS  tracking 

*  Convergence  of  the  iterative  corrections  of  the 
gravity  field  and  the  orbit  ? 

*  Simultaneous  estimation  of  the  orbit  and  the  grav¬ 
ity  field  ? 
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The  Use  of  6radiometers  in  Space  to  Monitor  Changes  in  the 
Earth's  Gravity  Field. 


Oscar  L.  Colombo.  University  of  Maryland  Astronomy  Program, 
Code  626,  NASA  Goddard  Space  Flight  Center,  Greenbelt,  Maryland 
20771. 


Tides,  and  a  variety  of  processes  of  non-tidal  nature  associated 
with  the  oceans,  the  cryosphere,  and  the  atmosphere,  exert  variable 
loads  on  the  solid  earth,  resulting  in  fluctuations  of  the  external 
gravitational  field.  A  gravity  gradiometer  in  orbit  can,  in  principle, 
monitor  those  changes  to  study  both  the  loading  phenomena  and  the 
mechanical  properties  of  the  earth’s  Interior  governing  the  response 
to  the  loading.  Current  space  techniques,  involving  laser  ranging  to 
spacecraft,  can  reveal  only  broad  zonal  features.  An  orbiting 
gradiometer  may  provide  a  more  complete  picture.  Given  sufficient 
accuracy,  and  enough  observing  time,  such  an  instrument  could 
reveal  the  geographical  distribution,  in  both  latitude  and  longitude, 
of  changes  that  occur  at  frequencies  ranging  from  daily  to  secular. 
The  gradients  of  such  gravitational  changes  have  most  of  their 
power  in  the  band  from  once  per  orbital  revolution  (100  minutes)  to 
once  per  tenth  of  revolution.  Because  of  their  long  wavelengths,  they 
can  be  sensed  at  much  higher  altitudes  than  the  sharper  signals  of 
crustal  origin  that  are  the  main  concern  of  missions  such  as  GRM  or 
Aristoteles.  Surface  forces  like  drag  are  much  weaker  and  less  of  a 
problem,  and  a  mission  may  last  for  several  years,  instead  of 
several  months.  Typically,  the  signals  are  of  the  order  of  10-7  E,  and 
Paik's  cryogenic  instrument  could  allow  their  resolution  at  the  1 
percent  level  after  one  year  of  continuous  observation.  The  bandwith 
of  the  GRM  device  (1  Hz)  is  much  larger  than  needed  for  this 
application.  However,  useful  life  can  be  seriously  limited  by  the 
gradual  boliling  off  of  the  liquid  He  coolant.  Perhaps  instruments  of 
a  different  kind,  able  to  operate  in  space  for  many  years,  may  be 
constructed  specially  for  sensing  the  long-wave  changes  in  gravity. 


THE  USE  OF  GRADIOMETERS  IN 
SPACE  FOR  MONITORING 
CHANGES  IN  THE  GRAVITY  FIELD 
OF  THE  EARTH 


UNIVERSITY  OF  MARYLAND  ASTRONOMY  PROGRAM 
(NASA  GODDARD  SFC.  CODE  626. 
GREENBELT.  MD.  20771.) 
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Wagnfji  and  McAdoo:  Time  Vakiations  in  the  Earto's  Gravity  Fifld 


«J79 


Vertical  Motion  from  Clacial  Rebound  (cm  a'^) 


Fig  5.  Pr»cnl-day  rale  of  change  in  Ihc  radial  position  uf  ihc  :>urfjcc  of  ihe  i>ulid  c.iriii  c  ilcul.iled  fruni  ^iniplilied 
model  of  poilgla'  d  rebound.  Unilt  are  centimeterii  per  year  Measured  uplifl  uilh  le^peii  lo  ^ea  level  [H’n  uiul  /Vlinr. 
I983J  ai  A,  Angermunn  River  (0.9  cm/yr.l;  R,  Richmond  Gulf  (1.25  cm, 'yr ).  C.  Churchill  II  Oun  'vr) 
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Table  1.  Some  Recent  Estimates  of  Temporal  Variations  in  Zonal  Harmonics  of  Uie  Earth’s 
Gravitational  Field 


Source 

Reference 

Variation 

AC20 

AC30 

•  earthquakes 

Chao  &  Gross 
(1987) 

nonperiodic 
w/long  period 
trend 

±5xl0*l3/yr 

±2xl0-13/yr 

•  deglaciation  rebound 
of  crust 

Yoder  etal. 
(1983)  ■ 

secular 

(observed) 

-3.0xl0*^Vyr 

n.a. 

Rubincam 

(1984) 

on  LAGEOS 

-2.6xlO^Vyr 

n.a. 

•  snow  cover 

Chao  et  al. 

periodic:  annual 

IxlO'^^amp 

6x10-11  amp 

(1987) 

semiannual 

3x10"^^  amp 

1x10-11  anip 

•  continental  drift 

Sconzo  (1980) 

secular 

±2xl0-^4/yr 

n.a. 

Greenland  moving  at 
lOcm/yr  in  latitude 
with  a  depth  of  im¬ 
mersion  of  50  km 

•  tidal  breaking 

Paddack  (1967)  secular 

<-5x10-13 

n.a. 

•  earth,  ocean  tides 

Christodoulidis  periodic 

-  variable  - 

-  variable  - 

et  al.  (1988) 
and  others 

(observed) 

nontidal  contributions  lumped  in 
tidal  recoveries  at  forcing 
frequencies 

•  air  pressure  & 

Gutierrez  & 

periodic:  annual 

IxlO'^amp 

n.a. 

groundwater 

Wilson  (1988) 

semiannual 

(shows  atmosphere/oceans  to  be 
~10  times  water  storage  at  annual 
and  ~3  times  at  semiannutd 
periods) 

1.5x10-19  amp  n.a. 

•  changes  in  sea 

Peltier  (1988) 

secular 

2xl0-li/yr 

n.a. 

due  to  ice  cap/ 
glacial  melting 

Yuen  et  tU. 
(1987) 

2  to  8xl0-l2/yr 

2to7xl0-l2/yr 

•  growth  of  the 
Antarctic  ice  sheet 

Yuen  et  al. 
(1987) 

secular 

5  to  lOxlO-12/yr 

6  to  1  lxlO-12/yr 

equivalent  to  drop 
in  sea  level  of 

0.3  mm/year 

•  continental  water 

Chao  (1988) 

periodic:  annual 

1.5x10-19  amp 

1.4x10*19  amp 

storage,  aquifers. 

semiannual 

5x10*11  amp 

4x10*11  amp 

lakes 

secular 

IxlO'l^/yr 

n.a. 
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POSSIBLE  TECHNIQUES  FOR  MEASURING  GRAVITY  CHANGES 
FROM  SPACE: 


SATELLITE  LASER  TRACKING 
PRINCIPLE:  maps  long-wave  zonal  signals  by 

DETECTING  LARGE  RESONANT  ORBITAL  PERTURBATIONS 

REQUIREMENTS:  one  lageos/starlette-type 

SATELLITE  FOR  EACH  TWO  ZONALS  (APPROX  .4  SATELLITES 
ROUGHLY  EQUISPACED  IN  INCLINATION  TO  RESOLVE  ZONAL 
CHANGES  TO  DEGREE  8.  OR  25  DEGREES  RESOLUTION) 

LIMITATIONS:  drag  corrupts  signals,  only  zonals 

RESOLVABLE.  SEVERAL  SPACECRAFT  NEEDED. 


6P5  TRACKING 

PR  I  NCI  PLE:  measures  whole  field  (zonal  and  non- 

ZONAL)  BY  TRACKING  OF  A  LOW  SPACECRAFT  CARRYING  A 
GPS  RECEIVER.  BY  SIGNALS  FROM  THE  GPS  SATELLITES. 

REQUIREMENTS:  additional  receivers  round  the 

WORLD  TO  CORRECT  CLOCK  ERRORS  IN  TRANSMITTERS  AND 
ORBITING  RECEIVER 

LIMITATIONS:  high  physical  stability  of  various 

COMPONENTS  NEEDED.  BUT  NO  GOOD  CONTROL  ON  CONDITIONS 
ON  GPS  SATELLITES.  OR  ON  GROUND  RECEIVERS;  DRAG. 
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6RM-TYPE  S ATELL I TE  S  AT  E  L  L  I  T  E 
TRACKING 


PRINCI  PLE:  MEASURES  WHOLE  FIELD  WITH  TWO  DRAG- 
FREE  SPACECRAFT  A  FEW  HUNDREDS  OF  KM  APPART  ON  SAME 
ORBIT.  TRACKING  EACH  OTHER  BY  TWO  WAY  DOPPLER/LASER. 

REQUIREMENTS:  drag  free  spacecraft,  high 

PHYSICAL  STABILITY  OF  COMPONENTS. 

LIMITATIONS:  requires  very  good  non- 

GRAVITATIONAL  FORCE  COMPENSATION  (DRAG.  RADIATION 
PRESSURE.  ETC.). 


CRYOGENIC  GRAVITY  6RADIQMETER 


PRINCIPLE:  measures  whole  field  by  sensing 

DIFFERENCE  MODE  BETWEEN  ALIGNED  ACCELEROMETERS  BY 
SENSING  WITH  S.Q.U.I.D.S.  THE  MAGNETIC  FLUX  DISPLACED  BY 
SUPERCONDUCTING  PROOF  MASSES. 

REQUIREMENTS:  similar  to  grm-type  sat.-sat. 

TRACKING. 

HIGH  PHYSICAL/MECHANICAL  STABILITY.  REJECTION  OF 
COMMON  MODE. 

LIMITATIONS:  self-gravitation.  vibrations. 

RESIDUAL  COMMON  MODE  ACCELERATIONS,  SCALE  FACTOR 
CALIBRATION. 
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CHARACTERISTICS  OF  A 
GRADIOMETER  MISSION  FOR 
MAPPING  TEMPORAL  CHANGES  IN 
GRAVITY 

ACCURACY:  10-5  TO  10-6  E  FOR  A  BANDWITH 
OF  0.01.  Hz 

BECAUSE  CHANGES  CANNOT  BE  MEASURED  DIRECTLY  ON 
EARTH.  A  GOOD  DEAL  OF  NEW  SCIENCE  CAN  BE  OBTAINED. 

SIGNALS  HAVE  LONG  SPATIAL  WAVE  LENGTHS.  SO  THERE  IS 
SLOW  ATTENUATION  WITH  ALTITUDE:  A  HIGH  ORBIT  (600- 
1000  KM)  CAN  BE  CHOSEN. 

WITH  ORBIT  600-1000  KM  HIGH:  MUCH  LESS  DRAG  THAN 
FOR  GRM  MISSION.  DRAG  FREE  SYSTEM  CAN  USE  HELIUM 
BOILOFF  OF  CRYOGENIC  GRADIOMETER  FOR  PROPULSION.  SO 
MUCH  LESS  WEIGHT  THAN  USING  HYDRAZINE  AT  200  KM  (MORE 
THAN  ONE  ORDER  OF  MAGNITUDE  LESS  FOR  PROPELLANT 
ALONE) 

A  SPACECRAFT  ALREADY  IN  DEVELOPMENT  (GP-B)  COULD  BE 
USED  (DRAG  FREE  USING  HELIUM  BOILOFF.  CRYOGENIC 
PAYLOAD.  ONE  AXIS  SPIN  STABILIZED.  PRECISE  ATTITUDE  IN 
INERTIAL  SPACE  DETERMINED  BY  TELESCOPES.  SPIN  MAY 
HELP  SEPARATE  SELF-GRAVITATION  AND  OTHER  SPURIOUS 
SIGNALS  FROM  THE  DESIRED  GRAVITATIONAL  INFORMATION. 

BECAUSE  OF  GREAT  SENSITIVITY  REQUIRED.  PROBLEMS  LIKE 
SELF-GRAVITATION  CAN  BE  DIFFICULT  TO  SOLVE. 


195 


Aj2  (CM  OF  WATER) 
1  CM 


A  ORBIT  (RESONANT,  IN  METERS) 


A  GRAVITY 
10**^®  G  1.5  M 


45335  45492  45649  45806  45963  46120  46277  46434 


Modified  Juliut  Date 


Fig.  1.  History  of  Starlcttc  node  residual  obtained  from  the 
three-year  continuous  orbit  using  the  nominal  force  and 


Tirtal  Perturbations  on  a  Satellite  Orbit 
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MEAN  SEA  SURFACE  TOPOGRAPHY 
(LEVITUS) 
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Figure  4-9.  SGG  Experiment  Module  concept  -  modified  GRM. 


Inversion  of  Airborne  Gravity  Gradient  Data, 

South-western  Oklahoma 

D.  W.  Vasco  (Center  for  Computatinal  Seismology  Lawrence 
Berkeley  Laboratory,  Department  of  Geology  and 
Geophysics,  University  of  California,  Berkeley,  CA 
94720;  415  486-7312) 

C.  L.  Taylor  (Geophysics  Laboratory,  Hanscom  AFB,  MA, 
01731:617  377-3078) 

We  present  a  preliminary  interperation  of  gravity 
gradient  anomalies.  The  diagonal  elements  of  the  gradient 
tensor,  as  recorded  by  the  Bell  airborne  Gravity  Gradient 
Survey  System  (GGSS),  are  used  to  compute  the  basement 
topography  in  south-western  Oklahoma.  This  is 
accomplished  through  a  non-linear  inverse  procedure 
based  on  the  conjugate  gradient  algorithm.  In  general  the 
resulting  model  contains  a  ridge  of  shallow  basement 
material  (<3.0  km)  trending  east  south-east.  This  ridge  is 
bounded  on  the  north  and  the  south  by  troughs  in  the 
basement  which  extend  as  deep  as  10.0  km.  The  gradient 
field  which  results  from  this  model  fits  most  of  the  GGSS 
observations  within  their  estimated  errors  of  12.0  E.  The 
depths  also  agree  with  a  set  of  available  oil  well  depths  to 
the  basement  and  with  inferred  faults  in  these  igneous 
rocks.  In  order  to  assess  the  derived  solution,  the  problem 
was  linearized  about  the  final  solution  and  linear 
parameter  resolution  and  parameter  covariances  were 
computed.  For  the  most  part  these  basement  depths  are 
well  resolved  and  the  resolution  matrix  is  diagonally 
dominant.  Futhermore,  the  parameter  standard  errors  are 
small,  the  majority  are  less  than  1.0  km.  Only  26 
parameters  out  of  98  have  errors  larger  than  this. 
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INVERSION  OF  AIRBORNE  GRAVITY 
GRADIENT  DATA,  SOUTH-WESTERN 
OKLAHOMA. 

D.  W.  VASCO  &  C.  L  TAYLOR 
GEOPHYSICS  LABORATORY  (AFSC) 
HANSCOM  AFB,  MA  01 731  -5000 
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For  a  prism  defined  by; 


a  <5  <a  ,|3  <P  ,Y  <5  <Y 

1  1  2  1  2  2  1  3  2 

where  ^  are  the  source  coordinates,  then  the 
diagonal  elements  of  the  gradient  tensor 
are: 


T  =pGarc1an(U/£r)| 

XX  •7  1 


T 

yy 


=  pGarctan(^  £  /£r)|  | 

13  2a 


T  pG  arctan(E  ^  r)| 


a 


where :  G  =  gravitational  constant 
p  =  density  contrast 
r  =  sourc  e-recei  ver  distan  ce 


The  objective  functional  is: 
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EAST (KM) 


NORTH 


0.00  3.50  7.00  10.50 


H.OO 


Condustans: 

*The  signal,  with  a  maximunn  of  57.5  E  is 
above  the  estimated  noise  level  of  1 2.0  E 
and  coherent  between  the  seven  tracks. 

*The  basement  model  presents  a  coherent 
structural  feature  trending  west-northwest  to 
east-southeast.  This  ridge  of  higher  density 
material  agrees  with  known  basement  faults. 

••• 

*The  model  resolution  of  a  majority  of  the 
prisms  is  adequate  and  the  standard  errors 
are  quite  low,  most  less  than  1 .0  Km. 

i 

I 

i 

!  I 

I 

I  *The  derived  solution  agrees  with  known 

!  basement  structure  and  available  oil  well 
data. 
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DEVELOPMENT  OF  A  MOBILE  GRAVITY  GRADIOMETER  FOR  GEOPHYSICAL  EXPLORATION 


FJ  van  Kann,  MJ  Buckingham,  MH  Dransfield,  AG  Mann, 
PJ  Turner,  R  Matthews,  RD  Penny  and  C  Edwards. 

Physics  Department, 

The  University  of  Western  Australia, 
Nedlands,  Western  Australia  6009. 


We  present  a  description  of  a  superconducting  gravity  gradiometer  designed  for 
geophysical  use.  The  initial  target  sensitivity  is  1  Eo/VHz  in  a  frequency  band  below 
1Hz. 

• 

The  OQR  instrument,  which  measures  an  off-diagonal  component  [xy]  of  the  gradient 
tensor,  consists  of  two  perpendicular  sensors  with  parallel  pivot  axes  aligned  along  the 
vertical  z-axis.  This  configuration  of  dual  Orthogonal  Quadrupole  Responders  enables 
rejection  of  angular  accelerations  about  the  pivot  axis.  Rotation  about  each  of  the  other 
two  axes  is  controlled  independently. 

Each  of  the  quadrupole  sensors  is  carefully  balanced  mechanically  at  room  temperature 
and,  since  the  pivot  is  integral  to  the  sensor,  this  balance  is  preserved  at  low 
temperature.  Residual  off-balance  compensation  and  matching  between  each  sensor  in  the 
pair  is  achieved  magnetically  using  superconducting  trim  coils. 

Motions  are  sensed  by  pairs  of  superconducting  pancake  coils  arranged  in  current 
differencing  configurations  with  SQUID  readouts.  Apart  from  that  desired,  the  signal 
from  the  primary  pair  of  coils  contains  small  residual  terms  resulting  from  the  common 
mode  accelerations  perpendicular  to  the  pivot  ax!.'>.  Signals  from  a  set  of  secondary  coils 
are  combined  passively  to  eliminate  the  effects  of  tne  common  mode  acceleration  vector. 
Residual  sensitivity  to  z-axis  angular  acceleration  is  treated  in  a  similar  way. 

Rotation  about  each  of  the  other  two  axes  is  measured  optically  and  controlled  by  a  servo 
referenced  to  a  room  temperature  Inertial  system.  The  latter  is  a  gimbalied  piatform 
stabilised  by  a  pair  of  phase  modulated  fibre  optic  gyros  to  about  2.10'^  (rad/sec)/VHz. 
The  cold  gradiometer  package  is  also  mounted  on  gimbals,  in  this  case  driven  by 
diamagnetic  actuators.  The  thermal  environment  in  which  the  package  is  mounted  is 
0  carefully  controlled  and  maintains  an  operating  temperature  constant  to  within  a  few  tens 

of  pK  at  about  5K  . 

Many  of  the  gradiometer's  features  have  been  proved  under  laboratory  conditions  and  we 
•  are  presently  engaged  in  testing  the  complete  package  prior  to  transferring  it  and  its 

support  systems  into  a  mobile  laboratory  in  readiness  for  field  trials.  Moving  base  tests 
are  scheduled  to  begin  before  the  end  of  the  year  at  the  Dongara  natural  gas  fields  3U0  km 
north  of  Perth  in  Western  Australia. 
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spring  balance  sensors 

(in  line  or  shear  gradients) 


same  mechanism  for  common  force  and  differential  force  stiffness 
=>  trade  off: 

low  stiffness  for  high  sensitivity 

high  stiffness  for  high  CMRR  (dynamic  range  problem) 


Beam  balance  sensors 

(shear  gradients) 


different  mechanisms  for  common  force  and  differential  force  stiffness 

•  low  differential  force  stiffness  »  1  Hz  high  sensitivity 

•  high  common  force  stifftiess  «  1  kHz  =>  high  CMRR  possible 

•  tune  CMRR  during  assembly  to  >  125  oB 
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t  =  %( gk(0)  Mj  +  rici(0)  Mji+ ....) 


gi  ^  gi 


Fij  =>  Fij  +  Rij  = 

Gij 

i.e.  symmetric  +  antisymmetric 

^  ^yy  "Mxx 

0 

°  1 

^  Fyx  "COyCOx 

T  = 
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Mxx  ■  Myy  ) 
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<  Cx  > 

0 

Mxx  +  Mxx 

0 

• 

«y 

.  0 

0 

Mxx  +  Myy  j 

,  Ox  > 

ti=^(12-b2)Gxy  m  =  mass,  1  *  length,  a  **  b  =  width 
Two  bars:- 

o«c(XA-tB)+  A(ta  +  xb)+  AAx*x+Myy+  (ASxx  -  ASyy)x  y 

I  I  I  I 

signal  common  mode  residual  dipole  induced  dipole 

TA  -  TB  =  2  ^  (P  -  b2)  (Fxy  -(OxWy) 

TA  +  TB  =  2^(12  +  b2)(a^) 

Q  =  earth  rotation  rate  =  7.3  10’^  rad/sec  =  IS^yhr 
5cOy,  S(0x  =  platfcim  angular  velocities 
cOxCOy  =>  Q  (sin0  cos(l).5a)y  -  cos9.5(Ox) 
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PIVOT  DEFORMATION 
UNDER  LOAD 


TORSION 

a 

L 

3  Hz 


STRETCH 

a  b 
L 

5  kHz 


s-BEND 

a  b3 
L3 

IkHz 


BAR  DEFORMATIONS  IN  A 
FORCE  FIELD  (IDEAL  PIVOT) 
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Angle  (deg) 
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Acceleration  amplitude  at  45' 


INFORMATION  FLOW 


Basic  detector  noise 


Equivalciu 

gradient 

noise 

(EO/VHz) 


Equivalent  gradient  noise  from  amplifier  l/f  and  Biownian  tenos 


Coloured  noise 

•  down  conversion 

•  thermal 

•  flux  creep 
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ROTATIONAL  STABILISATION 


SQUID 

ELECTRONICS 
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RESISTIVE  ANODE  DETECTOR 


ANGLE  SENSOR 


MOUNTED  ON  INERTIAL 
STABILISED  PLATFORM 


NOISE  LEVEL  O.l^rad/VHz  (6^deg/VHz) 
lOniliz  to  IkHz 

DYNAMIC  RANGE  O.Olrad  (0.6deg) 


COLD  MIRROR  ON  GRADIOMETER 


Optic  Gyro  -  StabiBzed  Platform 
Performance  Data 
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10  iniiuitON 


Set  up  Procedure 


XY  Gravity  Gradient: 

Spectral  Density  for  1  gm  cnf^  Density  Contrast 
and  UWA  Instrument  Sensitivity 


(zH/^/QH)  i-pods 
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DEVELOPMENT  OF  THE  MODEL  III 
SUPERCONDUCTING  GRAVITY  GRADIOMETER 


M.  V.  Moody,  Q.  Kong  and  H.  J.  Paik 

Department  of  Physics  and  Astronomy 
University  of  Maryland,  College  Park,  MD  20742 


The  development  of  a  three-axis  superconducting  gravity  gradiometer, 
SGG,  is  continuing  at  the  University  of  Maryland.  The  instrument  is  being 
developed  under  a  NASA  contract  for  the  purpose  of  precision  gravity  ex¬ 
periments  and  gravity  field  mapping  from  an  orbiting  platform.  Testing  of 
the  Model  III  SGG  has  recently  begun.  This  device  was  designed  to  meet 
the  sensitivity  requirements  of  NASA  for  a  global  gravity  mapping  mission 
(3  X  10-*  E 

The  SGG  utilizes  three  pairs  of  spring  mass  systems  in  which  proof  mass 
motion,  induced  by  a  gravitational  force  or  an  acceleration,  modulates  su¬ 
percurrents.  The  superconducting  circuits  are  configured  such  that  these 
supercurrents  are  passively  summed  and  differenced  before  being  measured 
by  SQUID  amplifiers.  Also,  in  order  to  operate  in  both  terrestrial  and  space 
environments,  the  proof  masses  in  the  SGG  use  a  superconducting  levitation 
scheme  which  has  minimal  effect  on  the  differential  mode  spring  constant. 

The  primary  enhancement  of  the  Model  III  over  previous  designs  is  the 
incorporation  of  a  passive  superconducting  negative  spring.  Using  the  nega¬ 
tive  spring  to  cancel  the  spring  constant  of  the  mechanical  spring,  the  noise 
contribution  of  the  SQUID  amplifier  can  be  suppressed.  The  results  of  the 
Model  III  SGG  tests  will  be  presented. 

Usbg  the  SGG  to  measure  the  Laplacian  of  the  gravitational  potential, 
a  composition  independent,  null  test  of  the  inverse  square  law  of  gravity  can 
be  performed.  Sensing  tilt  of  the  SGG  platform  with  a  laser  and  photodiode, 
we  have  demonstrated  that  tilt  is  the  primary  error  source  in  this  experiment 
when  using  a  1600  kg  pendulum  as  the  source.  Methods  for  reducing  this 
and  other  errors  using  a  three-axis  SGG  will  be  discussed. 


DEVELOPMENT  OF  THE 
MODEL  III 

SUPERCONDUCTING  GRAVITY 
GRADIOMETER 


M.  V.  Moody,  Q.  Kong  and  H.  J.  Paik 


Department  of  Physics  and  Astronomy 
University  of  Maryland,  College  Park,  MD  20742 
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ACCELEROMETER  CROSS  SECTION 
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SCHEMATIC  OF  MODEL  III  SGG 


CANTILEVER  SPRING 
-^-eZZ2m32Z  3-^  , 
r - y.  ./W'fffT'--!— , 


SENSING  CIRCUIT: 

Adjust  ratio  of  to  1^2  to  balance  out 
sensitivity  to  common-mode  accelerations. 
LEVTATiON  CIRCUIT: 

Energy  (0V2L)  is  constant  for  differential  motion, 
increases  only  common  mode 
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INTRINSIC  SPECTRAL  NOISE 


=  BROWNIAN  MOTION  +  AMPLIFIER 


FOR  BEST  COMMERCIALLY  AVAILABLE  SQUID: 
EaCF)  =  3  X  10'*“  J  Hz"' 


TO  REDUCE  AMPLIFIER  NOISE  CONTRIBUTION  LOWER  f,, 

1.  In  Qe  use  ’’push-pull”  levitation. 

f,  =  8  Hz.  S(f)  =  2  X  10‘®  E  Hz'* 

2,  Superconducting  negative  spring. 

f,  =  1  Hz.  S(f)  =  2  X  ID"*  E  Hz'* 
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SUPERCONDUCTING  NEGATIVE  SPRING 


aoons 

0.00«7 

0.00105 

0.00X0 

-I 

•ooowee 

aoo«tt 

0.001064 

aoot«3 

0.001082 

-X- 


ENERGY  =  Mon''Pv(0)/2V(x) 


NEGATIVE  SPRING 


NEGATIVE  +  LINEAR  SPRING 
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PRIMARY  ERROR  SOURCES 


=  r 

+  CENTRIFUGAL  ACCELERATION 
[1  “ 

+  COMMON-MODE  ACCELERATION  (TILT) 
(1//)((5n.^  +  h,A)-0(/)xgE 


+  ANGULAR  ACCELERATION 


CALIBRATION  AND  ERROR  COEFFICIENTS 


ADJUST  DRIVE  CURRENT  IN  TRANSDUCERS  TO  OBTAIN 
X  TILT.  Y  TILT  OR  VERTICAL  SHAKING. 

COMMON-MODE  CALIBRATION 

g  =  A-BxgE 

GRADIOMETER  CALIBRATION 

r(2f)  =  [1  -  (A*fi)W(f) 

MEASURE  MISALIGNMENT 

(5r(f)  =  -(1//)(5n.^‘exgE  +  ($n.^xh‘27Tf0 

(5n_^  =  3.4  X  10’^  (adjusted  at  room  temperature) 

(5n+^  =  5.0  X  10"^  (not  adjusted) 

DETERMINE  CENTRIFUGAL  ACCELERATION  IN  TWO  DIMENSIONS 
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NOISE  IN  THE  SENSING  CIRCUIT 


FLUX  LEAKAGE 
out:  0(i)  -  -  a(l  - 

in:  =  a0(i  - 

T  =  50  s,  a  =  2  X  10’®,  with  B  =  0.5  tesla 


BALANCE  BOTH  SENSING  CIRCUITS 
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Composition  Independent 

Null  Test  of  the  Inverse  Square  Law  of  Gravitation 


Test  for  non-Newtonian  potential  of  the  form 
$(r)  =  -G®(1  + 


TILT  IS  THE  PRIMARY  ERROR  SOURCE 
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REDUCING  SENSITIVITY  TO  TILT 
IN  A  SINGLE-AXIS  GRADIOMETER 

6r(cj)  =  6n_^-0(cj)xg£; 

=  0(cj)-gEX(5n_^ 

-  1^2  -L  ^ 

1.  ALIGN  f\  HORIZONTAL 

2.  ROTATE  ABOUT  A  TILL  ($n_^  II  g^ 
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Removal  of  tilt  and  scale  factor  mismatch  errors 
with  a  three  axis  gradiometer. 


r;,  =  ar,j+  h,n,-^x4 

r;3=^r33+  h3n3-i5x 

^33  =  7^33+  Se 


£C'-  r*tC''-fcC 


s  a 


B 


Rotate  gradiometer  120  degrees  twice  and  sum, 


'C  'C) 


DEVELOPMENT  OF  A  SUPERCONDUCTING  SIX-AXIS  ACCELEROMETER 


E.  R.  Canavan,  H.  J.  Paik,  and  J.  W.  Parke 

Department  of  Physics  and  Astronomy 
University  of  Maryland,  College  Park,  MD  20742 


The  three-axis  superconducting  gravity  gradiometer  being  developed  at 
Maryland  for  an  orbiting  gravity  mapper  requires  very  precise  platform 
stabilization,  particularly  against  angular  motion  noise.  The  key  component  of 
the  stabilized  platform  is  a  superconducting  six-axis  accelerometer.  The 
accelei’ometer  can  also  function  as  a  complete  inertial  navigation  system,  anrl 
with  the  gradiometer  it  forms  a  gradiometer-aided  inertial  navigation  system. 

The  device  senses  the  motion  of  a  single  levitated  niobium  proof  mass  with 
respect  to  its  housing  using  superconducting  AC  inductance  bridges  and  a 
SQUID  amplifier.  The  proof  mass,  composed  of  three  intersecting  square 
slabs,  fits  inside  a  housing  of  complementary  shape  formed  by  8  titanium  cubes 
mounted  in  the  corners  of  a  large  hollow  cube.  The  face  of  each  titanium  cube 
adjacent  to  the  proof  mass  holds  a  levitation  and  a  sensing  coil.  The  24 
levitation  and  the  24  sensing  coils  are  connected  to  form  circuits  that  provide 
levitation  and  sense  displacement  in  each  of  the  6  degrees  of  freedom. 

The  first  prototype  of  the  device  has  been  built  and  operated  successfully. 
The  measured  values  for  resonance  frequency,  sensitivity,  and  other 
parameters  matcfi  very  well  to  those  given  by  a  detailed  analytical  model.  The 
model  predicts  that  by  optimizing  electro-mechanical  coupling,  which  al 
present  is  small,  and  using  a  better  SQUID,  the  accelerometer  should  be  able  to 
achieve  a  base  noise  level  of  10'^^  g/VHz  and  rad/s^/VHz.  Larger  coupling 
should  be  achieved  in  a  prototype  under  development. 
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Development  of  a 
Superconducting 

Six  Axis 
Accelerometer 

E.R.  Canavan,  H.J.  Paik,  &  J.W.  Parke, 
University  of  Maryland, 

College  Park,  MD  20742 
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To  develop  an  accelerometer  that  is: 

•Extremely  sensitive 

•Compact 

•Measures  all  6  degrees  of  freedom 

•  Compatible  with  the  SGG 

Principle  of  Operation: 

•  Single  magnetically  levitated  mass 

=>  responds  in  all  degrees  of 
freedom: 

q  ->q  ,  q  =  {  r^,  ry,  r„  0^,  Gy,  0, } 

•  Displacements  alter  the  inductance 

of  24  coils  surrounding  the  mass: 

L  =  Lo  +  Aq  +  0(  q' ) 

•  Coils  are  arranged  into  6  AC 

inductance  bridges,  each  sensitive 
to  motion  in  a  different  degree  of 
freedom. 


Sensing  Circuit 


Oscillator 


Upper 

coils 


SQUID 


Lower  coils 


Operation: 


•Circuit  analysis  gives: 


(L2L3  -  L1L4) 

^OSC _ 

(Li  +  LzJLs  +  LiI  +  LsqXLi 


•By  geometry, 

Li,  L4  =  Lq  -  Astx  ;  Lg,  L3  =  Lq  +  Asr^ 


isQ  = 


/  A  • 

As  Ipsc 

^Lo  +  Lsq  , 


=»  SQUID  output,  after  demodulation, 
is  proportional  to  r^. 
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Levitation  Circuit 


Coils 


Proof  Mass 


from  Feedback  ^ 

controller  Heat  Switch 

Pulsing  heat  switch  while  applying 
current  II  traps  II  in  the  loop. 
Choose  II  to  minimize  q. 


V  = 


=  (*1-0  ^  Lfllf,  _  2^ ^  8gg^ 

2  4Lo  +  Lf 


=5fDC  =  2ALlL  k  = 
•Feedback  current  adds  to  II. 
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Multiplexed  Operation 

All  6  sensing  circuits  are  connected  in 
series  with  a  single  SQUID.  Each 
bridge  is  driven  at  a  different 
frequency  Oj  and  the  output  of  the 

SQUID  is  fed  to  6  lock-in 
amplifiers  where  the  signals  from 
the  6  bridges  are  demodulated. 
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Materials  Considerations 


•Coil  Forms 

Material:  Ti6A14V 

Problem:  sensitive  to  heat 

treatment;  in  our  case,Tg>4.2K. 

=>  Need  temperature  controller  to 
maintain  SSA  above  T^. 

•  Superconducting  wire 
Material:  NbTi 

Problem:  Alloy  superconductors  are 
Type-II  drift  in  II  due  to  flux 

creep  =>  low  frequency  noise. 

Originally  used  Type-I  Nb,  but  had 
significant  occurance  of  thermal 
stress  breakage.  (Improper 
drawing  process.) 
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ANALYriCAt-  Mo&Ct- 
OGENEf^M  EQUATIONiOF 
MOTION^ 

Pnoof  /IASS  MTIOII  w.  n.r. 

HOUiwi  PCS6RIBet>6Y  A'cei,/^' 

Hoojia/S  MoriM  ma.T.  /NEATiAL 

PAA/ie  PEStRiBep  BY  A"  r**. 

S>  V’  •  «>/'x(AV')*A1''**  r* 
io“  •  w"*  A'w' 


wdAJS'rKtAtYlAXA  TH-e  LA6RAAI6(A//  A»/D 

l>£AlVIA/6  the  of 

--r*-  r^-  w'x(<w»'V/‘^]-wV^ 
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rne  Equ/^t/oa/s 

hOrtOM  ARE  AJOAiLlA/EAA. 
TO  Fmr  ofiptR! 

**  A 

•  -H 
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2.CMCULATI0N  OF  pOfEA/nAL 

CMCtiLhJ-K  V(9,/A  As  M 
lurAobucTro/u^  Bur  iajcludc.  i***  0  RPe/t 
TpBAiSi  L»  L**Ax-|x‘-4  B' 
CAtCUCATi/vlC  ^  TMC  LEVlThTiOX 
AM^  S6AJSfAJ6  CfACUirs  EACf^ 
AXf5  V  S{J/iJ1IAj6  • 

y  s  Vo  •  f>t  ( /«  ♦  Ty  ♦  /^ ) 

*  ^  (/y-/;)) 

♦  iC'^i^hiKOrV;  -♦ 

«>  toHiAif^Ajr  sfMj^c  ccMiT/^A^rs: 
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3.  ChucuLt^re.  Tmajsw  FvMtT/OAi 

rj  t  0  «  A; 

♦  w^'  d,-  » ot;  wljf/ie,  ex/»/iai6^ 


•c,«  0y  (o* **^ 

^  H/^aIT  to  U&S  COAjT/lOt-tf/C 

C^LeOuATlM6  5EAiSIAJ6  Cl^cOir  TA^/JSFC/ 
I^UAjtTfOAJ  ^  (as  lA)  lAiTROtjOdr/ 0>u) 
CO/^  CiaJiaJC  vv/ith  RBOVC- 

<A. 


ui  S  s  I  <J\.«  I 

WOTe:HR‘/E  ftt^PED  Vei^i^nr-PEPSMPSAJT  ^A^pi/Jf 


H.  WlNinUn  DETECTABLE  AaClEAATfO/ 


TV^  F(;AiDA/1£A)TAL  NOIse 

•BaowimiAa^  nonoAi  NOtse: 

ACCEL6AATI0M  sftcUH.  x>en^nrj 

S|T8  2isl±2L 
*  ^  0( 

*  SQUID  A/tN.ineA  Noise: 


(/siN6  H^; ,  cA<a  ^ctArrc  iMfwr 
cuA^eAir  NOfse  srecTAAL  peA)«in^S| 

To  B^ut^HtAJT  kcceLEKhrm 

SfZCTAAL  P6A/SITV'|  * 

niNi/iv/i  ^ersc-TABLe 


wHeAE*  emcaoh  coQfLtAJo  coen^iciejjy 
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FOA  THfc  AflOULMA  P £0^665  df 


FACEDO/t^  O^TiklAii 


VMN&AC 
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COMPA'RISON  OF  HODEL  WITH 
EXTERItlENTAL  'RESULTS 

I)  i^esonant  frequences 

TO  CALCULATE  NEED  Al.S, 

COMPUTE  inductance 'PA'RAMETERS  FRO^ 
FORCE  ‘BETWEEN  SET  OF  CONCENTRIC 
LOOPS  ANT)  their  /HAQE  CURRENTS 


EXPAND  AND  NUME'R.'CALLY  INTE6RATE 

A  > 

SIM1LAT2LY  TORQUE  EQUATION 
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5UBSTITUTIN6  ?;4T0  EQUATION  FOR, 

K/r*  (SET 

U>,  .  '  KkiS 

‘  '  m  ^  ’  X 


TCALC 

fexp  L' 

Px 

21.  Z 

2J.G 

Py 

20. 

21. S 

Pz 

17.8 

10,2 

13  3 

15,01 

I3.<^ 

IS2[ 

i^Z 

m.5 

|6.bJ 
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MORE  REPENOENT 

ON  p 


2)  SENSITIVITY 


TO  CALCULATE  Haxi.MOSTKNOWli 

■PROBLEM  -  BRIDGE  ‘DRIVEN  WITH 

TANK  CIRCUIT  (TO  IMPROVE 
SAIN, -REDUCE  •RF  INTERFERENCE) 

SOLUTION  -  MEASURE  Z(w) 

-»  CALCULATE  CIRCUIT  SAIN 


COMPUTE  I5 

■RESULTS 


^Hqxi^cxp 

y 

3.1  X  iqf 

3.5-‘IO^' 

Y 

^  10, 

8.8«io; 

1 

l.€x|0* 

Z.U-IO  J 

Oy 

2.1  X 10^ 

2. 1  -•  10^ 

^Y 

60x10 

S.2»I0 

6z 

Uxlo' 

k.  1  *  10? 

REASONABLY  GOOD  FIT 
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a)  MINIMUM  SI&NAL 

FOR  QUANTUH  T)eS>6N  5®OI0S 

Es*IO*“3/Hi  sl"  =  10'“^.//^ 

ACTUALLY  SEE  TWIS  LEVEL  IN  MEASUREHENT 


POWER  SPECl 
1.0  — 


20Ayg _  0%0vlp  Hann 


50  .  O 


Fxd  ^  O  Hz  Y  Q 

EXPERIHENm  "P^  • 

CP^^HOO, 


mJWw\u 

200 


WkW 

y  3.2*10’’! 

Y  1.0  *I0”* 

Z  S.M  *  lO"’® 

0x  M.e^io;’ 

Oy  «.7x|0\ 


/cHqkO 


V « ««  ^Hooct 
M.QylO"**)  % 

6,7*  IO*/»2i 

9  I  S6C 


2.2*10 


^  -C 
3.b*lO  - 

2.0  *  loy 
2.0*  10’® 


Vh> 


OTHER  ■RESULTS 

I]  FRESUeNCY  NOISE 


FREQUENCY  ( Hz ) 


NOTE  •  “surf"  FEAK  CLEARLY  VISIBLE 
SEISHIC  (WISE  SHOUU)  NAVE  HINIHUh 

AT  io-"'inj,) 

^  Low  FReO.  NOISE  UNITS  SM 

BELOW  ~lor^H2 

LOVJ  FR^OI*  noise  CC  V 


o  in 


l)  FEED  -BACK 

NEED  FOR  CONTROLLER* 

•  LINEARIZE  OUTPUT 

•  REDUCE  CROSS  COUPLIN6 

•  INCREASE  DYNAMIC  MNGE 
USE  TID  CONTROLLER 


OPEN 

^‘LOOP 


CLOSED 

^LOOP 


FEEDBACK  WIDENS  DANOWIOTH 


CROSS  COUPLIN6 


STRONG  CROSS  COUPLINGs^MAY  NeeO  Hiho 

CONTROU0? 

NEAGORE  on  and  off  “DIAGONAL  RE^POMSE 


NOISE 

Yl-^OOl 


-360 


Noise-*' 

yo-*oo 


* 


FIND  OFF- DIAGONAL  GAIN  >10«IB  DELOW 


diagonal  ^system  -DIAGONALLT  T)0MINANT 

^CAN  use  SISO  CONTROLLER 


OPTIMIZATION 
•  ULTIMATE  performance- 

Prt  -j^rfKgT 

“0*  M  Vq,  *  P<r  J 


C  Pr)oTT  '  -i 

FOR  PEST  SQUIOj,  Eg'ZO'A 

*I.H  *‘10'^'!^^ 

Th? 

(Qt‘10*,  T»UK) 
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•IMPROVIMG  P  2 
<®»  -  »  /  IrAs  \  Lsq 


OSINS  SePCRATE  SENSING  CIRCUITS' 

GU  -^Ls 


IMPeOENCE  /lATCWNG  TRANSFORHeRS-' 


Lsq  “^Ls 


decrease  :  h,  Wpid 


INCREASE  -  A  ’♦INCREASE  TURNS 

DENSITY,  AREA 

Tr  ^IMPROVE  COIL 
MATCHING 

^single  uter  coils. 


MODEL  I  *  LARGE,  single  UYER 

SENSE/LEV.  Coil 


for  Ir'IAHP 
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SUPERCONDUCTING  GRAVITY  GRADIOMETER  MISSION  -  AN  OVERVIEW 


Ho  Jung  Paik 

Department  of  Physics  and  Astronomy 
University  of  Maryland,  College  Park,  Maryland  20742 


Two  dedicated  space  missions  proposed  for  the  1990’ s  hold  the  promise  of 
providing  data  for  recovering  the  Earth’s  gravity  anomaly  with  unprecedented 
accuracy  and  resolution:  the  Aristoles  Mission  and  the  Superconducting 
Gravity  Gradiometer  Mission  (SGGM).  SGGM,  the  more  ambitious  of  the  two,  aims 
at  recovering  the  global  gravity  field  to  a  precision  of  two  to  three  mgal 
with  a  resolution  of  50  km. 

The  instrument  package  of  SGGM  is  a  three-axis  gravity  gradiometer  which 
is  integrated  to  a  six-axis  accelerometer  for  active  platform  control.  The 
intrinsic  sensitivity  of  the  gradiometer  is  10*“*  E  and  that  of  the 

•13  *1/2  -11  -2 

accelerometer  is  10  g  Hz  in  linear  acceleration  and  10  rad  sec 
-1/2  ^ 

Hz  in  angular  acceleration.  'While  precise  attitude  control  of  the 
Experiment  Module  is  essential  to  mission  success  and  is  also  technically  most 
challenging,  pointing  accuracy  and  disturbance  isolation  requirements  of  SGGM 
are  less  stringent  compared  to  that  of  other  missions,  such  as  Hubble  Space 
Telescope  (HST)  and  Gravity  Probe-B  (GP-B).  Thus,  they  are  within  the  reach 
of  technologies  of  the  1990' s. 

In  the  recently  completed  Phase  A  study,  the  SGGM  study  team  addressed 
the  problem  of  scientific  requirements  and  mission  feasibility.  At  the 
University  of  Maryland,  prototypes  of  the  three-axis  gradiometer  and  the 
six-axis  accelerator  are  being  fabricated,  improved  and  tested.  The  actual 
mission  hopefully  will  take  place  before  the  year  2000. 


SUPERCONDUCTING  GRAVITY  GRADIOMETER  MISSION 
-  AN  OVERVIEW 


Ho  Jung  Paik 

Department  of  Physics  and  Astronomy 
University  of  Maryland,  College  Park,  MD  20742 


1.  SCIENCE  OBJECTIVES 

2.  SUPERCONDUCTING  GRAVITY  GRADIOMETER 

3.  SPACECRAFT  AND  ORBIT 

4.  DEVELOPMENT  SCHEDULE 

5.  CRYOGENIC  REQUIREMENTS 


October  13,  1989 

17th  Gravity  Gradiometer  Conference 
Hanscom  AFB,  MA  01731 
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SPECTRUM  OF  THE  VERTtCAL  GRAVITY 
GRADIENT  (EU  PER  COEFFICIENT) 


Horizontal  Resolution 


2.  S(4P€rt*nducHrty  ^rduf/cmef^f 


f)  Cf*r0f»0»  tit  LHe  t-em/KrtfuttS  (T^hT'^'iftm.K) 
4#m/  Bi^wnt^in  ntcffion  t^efsc 

^Hi^h  StAk!l;*y 

2)  "TmtstlucBtt  tint!  dt^9r*ti*i^^  (ty  ptrUittnt  «Mrrt^ 

^StAl*  factor  aaJ  itulf  stabifffy 
l-arff€  d/ c.  taaife 

i)  AnpHti tato0i  hy  SAUtD 
Hi^h  S'ens/t'Cif^ 

La/J<  dynamic  tanat 

H)  ^af'truinditC’Hin  nt^atra  tyoinj 
Ai^h  S'enSi^¥f^ 

S')  Sufcrflaid/namj  hel>'ttm>  baH^i 
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Intrinsic  Spectral  A^aise 

iif  (iwf,) 

Q(f)  ^  ifin 

Design  parameters 

prcaf  mass 

m 

uk, 

base  line 

i 

CJ9  m 

resonance  fre^uenc^ 

f. 

<lHt 

temperatu  re 

T 

UK 

puald^  factor 

a 

amplifier  net  sc 
(sue  die  ^^UID) 

r* 

3t/e"° 

enerc^u^  ceuplinc^  fatter 

fin 

Gravt  i(j  Gradient  ^oi$e 

Wiihoui  spring  Wtih  ne^atine  Sjon'n^(t,-l 

ff)  =  I  >cio-’ E ti  O.lHi  E  tfo.lHi 


Goal  se(  bij  I9S'3  watksha^  Sx/o'*  £  //i  ^ 


SIGNAL  FREQUENCY 


PARAMETEn  I  ERROR  MECHANISM  b  ORIENTATION  I  REQUIRED  CONTROL/KNOWLEOOE 


(a)  Proof  Mass  (Niobium) 


(c)  Assembly  Drawing 


(d)  Cross  Section  View 


Fig.  5.  Six-Axis  Superconducting  Accelerometer 


SIGNAL  FREQUENCY  f(H2) 
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SGGM  Spacecraft 
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